Search for: All records

ABSTRACT The multimessenger combination of gravitational waves (GWs) from merging massive black hole binaries (MBHBs) and the electromagnetic (EM) counterpart from the surrounding circumbinary disc (CBD) will open avenues to new scientific pursuits. In order to realize this science, we need to correctly localize the host galaxy of the merging MBHB. Multiwavelength, time-dependent EM signatures can greatly facilitate the identification of the unique EM counterpart among many sources in LISA’s localization volume. To this end, we studied merging unequal-mass MBHBs embedded in a CBD using high-resolution 2D simulations, with a $$\Gamma$$-law equation of state, incorporating viscous heating, shock heating, and radiative cooling. We simulate each binary starting from before it decouples from the CBD until just after the merger. We compute EM signatures and identify distinct features before, during, and after the merger. We corroborate previous findings of a several orders of magnitude drop in the thermal X-ray luminosity near the time of merger, but with delayed timing compared to an equal-mass system. The source remains X-ray dark for hours post-merger. Our main results are a potential new signature of a sharp spike in the thermal X-ray emission just before the tell-tale steep drop occurs. This feature may further help to identify EM counterparts of LISA’s unequal MBHBs before merger without the need for extensive pre-merger monitoring. Additionally, we find a role-reversal in which the primary out-accretes the secondary during late inspiral, which may diminish signatures originating from Doppler modulation. 

Abstract The physics of turbulence in magnetized plasmas remains an unresolved problem. The most poorly understood aspect is intermittency—spatiotemporal fluctuations superimposed on the self-similar turbulent motions. We employ a novel machine learning analysis technique to segment turbulent flow structures into distinct clusters based on statistical similarities across multiple physical features. We apply this technique to kinetic simulations of decaying (freely evolving) and driven (forced) turbulence in a strongly magnetized pair-plasma environment, and find that the previously identified intermittent fluctuations consist of two distinct clusters: (i) current sheets, thin slabs of electric current between merging flux ropes, and; (ii) double sheets, pairs of oppositely polarized current slabs, possibly generated by two nonlinearly interacting Alfvén-wave packets. The distinction is crucial for the construction of realistic turbulence subgrid models. 

Abstract The black hole at the center of M87 is observed to flare regularly in the very high-energy (VHE) band, with photon energies ≳100 GeV. The rapid variability, which can be as short as 2 days in the VHE lightcurve, constrains some of the flares to originate close to the black hole. Magnetic reconnection is a promising candidate for explaining the flares, where the VHE emission comes from background soft photons that inverse Compton scatter off of high-energy electron–positron pairs in the reconnecting current sheet. In this work, we ray trace photons from a current sheet near the black hole event horizon during a flux eruption in a magnetically arrested state in a general relativistic magnetohydrodynamics simulation. We incorporate beaming of the Compton up-scattered photons, based on results from radiative kinetic simulations of relativistic reconnection. We then construct VHE lightcurves that account for the dynamics of the current sheet and lensing from general-relativistic effects. We find that most of the flux originates in the inner 5 gravitational radii, and beaming is essential to explain the observed flux from the strongest VHE flares. The ray traced lightcurves show features resulting from the changing volume of the reconnecting current sheet on timescales that can be consistent with observations. Furthermore, we find that the amount of beaming depends strongly on two effects: (i) the current sheet inclination with respect to the observer and (ii) the anisotropy in the direction of motion of the accelerated particles. 

Abstract Using two-dimensional general relativistic resistive magnetohydrodynamic simulations, we investigate the properties of the sheath separating the black hole jet from the surrounding medium. We find that the electromagnetic power flowing through the jet sheath is comparable to the overall accretion power of the black hole. The sheath is an important site of energy dissipation as revealed by the copious appearance of reconnection layers and plasmoid chains. About 20% of the sheath power is dissipated between 2 and 10 gravitational radii. The plasma in the dissipative sheath moves along a nearly paraboloidal surface with transrelativistic bulk motions dominated by the radial component, whose dimensionless 4-velocity is ∼1.2 ± 0.5. In the frame moving with the mean (radially dependent) velocity, the distribution of stochastic bulk motions resembles a Maxwellian with an “effective bulk temperature” of ∼100 keV. Scaling the global simulation to Cygnus X-1 parameters gives a rough estimate of the Thomson optical depth across the jet sheath, ∼0.01–0.1, and it may increase in future magnetohydrodynamic simulations with self-consistent radiative losses. These properties suggest that the dissipative jet sheath may be a viable “coronal” region, capable of upscattering seed soft photons into a hard, nonthermal tail, as seen during the hard states of X-ray binaries and active galactic nuclei. 

Abstract We study the linear stability of a planar interface separating two fluids in relative motion, focusing on conditions appropriate for the boundaries of relativistic jets. The jet is magnetically dominated, whereas the ambient wind is gas-pressure-dominated. We derive the most general form of the dispersion relation and provide an analytical approximation of its solution for an ambient sound speed much smaller than the jet Alfvén speedv_A, as appropriate for realistic systems. The stability properties are chiefly determined by the angleψbetween the wavevector and the jet magnetic field. Forψ=π/2, magnetic tension plays no role, and our solution resembles the one of a gas-pressure-dominated jet. Here, only sub-Alfvénic jets are unstable ( $0<M_e\equiv \left(v/v_{\mathrm{A}}\right)\cos \theta <1$, wherevis the shear velocity andθthe angle between the velocity and the wavevector). Forψ= 0, the free energy in the velocity shear needs to overcome the magnetic tension, and only super-Alfvénic jets are unstable ( $1<M_e<\sqrt{{(1+{\mathrm{\Gamma }}_w^2)/[1+(v_{\mathrm{A}}/c)}^2{\mathrm{\Gamma }}_w^2]}$, with Γ_wthe wind adiabatic index). Our results have important implications for the propagation and emission of relativistic magnetized jets. 

ABSTRACT The upcoming Laser Interferometer Space Antenna (LISA) is expected to detect gravitational waves (GWs) from massive black hole binaries (MBHB). Finding the electromagnetic (EM) counterparts for these GW events will be crucial for understanding how and where MBHBs merge, measuring their redshifts, constraining the Hubble constant and the graviton mass, and for other novel science applications. However, due to poor GW sky localization, multiwavelength, time-dependent EM models are needed to identify the right host galaxy. We studied merging MBHBs embedded in a circumbinary disc (CBD) using high-resolution two-dimensional simulations, with a Γ-law equation of state, incorporating viscous heating, shock heating, and radiative cooling. We simulate the binary from large separation until after merger, allowing us to model the decoupling of the binary from the CBD. We compute the EM signatures and identify distinct features before, during, and after the merger. Our main result is a multiband EM signature: we find that the MBHB produces strong thermal X-ray emission until 1–2 d prior to the merger. However, as the binary decouples from the CBD, the X-ray-bright minidiscs rapidly shrink in size, become disrupted, and the accretion rate drops precipitously. As a result, the thermal X-ray luminosity drops by orders of magnitude, and the source remains X-ray dark for several days, regardless of any post-merger effects such as GW recoil or mass-loss. Looking for the abrupt spectral change where the thermal X-ray disappears is a tell-tale EM signature of LISA mergers that does not require extensive pre-merger monitoring. 

ABSTRACT We study the linear stability of a planar interface separating two fluids in relative motion, focusing on the symmetric configuration where the two fluids have the same properties (density, temperature, magnetic field strength, and direction). We consider the most general case with arbitrary sound speed cs, Alfvén speed vA, and magnetic field orientation. For the instability associated with the fast mode, we find that the lower bound of unstable shear velocities is set by the requirement that the projection of the velocity on to the fluid-frame wavevector is larger than the projection of the Alfvén speed on to the same direction, i.e. shear should overcome the effect of magnetic tension. In the frame where the two fluids move in opposite directions with equal speed v, the upper bound of unstable velocities corresponds to an effective relativistic Mach number $$M_{\rm re}\equiv v/v_{\rm {f}\perp }\sqrt{(1-v_{\rm {f}\perp }^2)/(1-v^2)} \cos \theta =\sqrt{2}$$, where $$v_{\rm {f}\perp }=[v_{\rm {A}}^2+c_{\rm s}^2(1-v_{\rm {A}}^2)]^{1/2}$$ is the fast speed assuming a magnetic field perpendicular to the wavevector (here, all velocities are in units of the speed of light), and θ is the laboratory-frame angle between the flow velocity and the wavevector projection on to the shear interface. Our results have implications for shear flows in the magnetospheres of neutron stars and black holes – both for single objects and for merging binaries – where the Alfvén speed may approach the speed of light.